Phonon modes at finite temperature in graphene and single-walled carbon nanotubes 
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The phonon modes at finite temperature in graphene and single-walled carbon nanotubes 
(SWCNT) are investigated by the mode projection technique combined with molecular dynam- 
ics simulation. It is found that the quadratic phonon spectrum of the flexural mode in graphene is 
renormalized into linear at finite temperature by the phonon-phonon scattering. Possible influence 
of this renormalization on the electron-phonon interaction and phonon thermal transport is ana- 
lyzed. For the SWCNT, however the quadratic property of the phonon spectrum for the flexural 
mode can survive at finite temperature under protection of the two-fold degeneracy of the flexural 
mode in SWCNT due to its quasi-one-dimensional structure. The frequency and life time of op- 
tical phonon modes are also studied at different temperatures. In particular, the life time of the 
in-plane optical phonons in graphene and the axial optical phonon in SWCNT can be described by 
an uniform formula r(T) = 560/T, which coincides with the experimental results. 

PACS numbers: 63.22.Rc, 63.20.kg, 02.70.Ns, 62.25.Jk 
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I. INTRODUCTION 

The phonon spectrum in graphene and single-walled 
carbon nanotubes (SWCNT) is a fundamental and im- 
portant issue. Phonons can considerably affect elec- 
tron conduction and thermal transport in carbon-based 
nano-devices, and the electron-phonon interaction in- 
duced Joule heating is one of the most common origins 
for the breakdown of electronic nano-devices. Consid- 
ering the importance of phonon modes in graphene and 
SWCNT, there has been lots of experimental and theo- 
retical works on this topic (for review, see, e.g. Ref.[l|). In 
the theoretical models satisfying the rigid rotational sym- 
metry, a special phonon mode with quadratic spectrum 
can be found which is the so-called flexural mode^— 
In graphene, flexural mode is the out-of-plane acoustic 
phonon mode. It is non-degenerate due to the sheet 
configuration of the graphene. The flexural mode of 
SWCNT is a two-fold degenerate transverse acoustic 
phonon mode. This two-fold degeneracy is due to the 
quasi-one-dimensional structure of the SWCNT where 
the two lateral directions are equivalent to each other. 
The quadratic property of the phonon spectrum ensures 
the flexural mode to be the lowest-frequency phonon in 
the long wave length limit, thus it will be the first mode 
to be thermally excited. As a result of its quadratic 
phonon spectrum, the flexural mode plays a special role 
in the electron-phonon interaction. In Ref. |6|, Mariani 
et at found that the electron scattering rate from flexural 
phonon modes will diverge logarithmically with temper- 
ature due to the quadratic phonon spectrum of flexural 
mode. In the ballistic thermal transport in graphene, the 
transmission function for flexural mode is T[oj] oc s/lj due 
to the quadratic phonon spectrumi^r— As a result, the 
thermal conductance in low temperature region is domi- 
naed by the flexural mode and depends on temperature 
as k oc T 1 - 5 . 

Since the qudratic property of the phonon spectrum for 



flexural mode is important and experiments are carried 
out at finite temperatures, it is necessary to investigate 
carefully whether the quadratic property is preserved at 
finite temperature or not. At finite temperature, phonons 
can have interaction with electrons, phonons, point de- 
fects, and etc. In present paper, we focus on the effect 
of the phonon-phonon scattering (PPS) on the phonon 
spectrum. The non-equilibrium Green's function is a 
standard yet complicate method to study properties of 
the phonon modes with PPS at finite temperature.— In 
another approach, the PPS can be included in empirical 
potentials, and its effect on the phonon spectrum can be 
extracted by molecular dynamics (MD) simulation and 
velocity correlation function.— ~— Using this method in 
SWCNT, it was found that the flexural mode still has 
a quadratic phonon spectrum*^ It will be interesting to 
investigate the PPS effect on the phonon spectrum in 
graphene by running MD simulation, since the structures 
of graphene and SWCNT are essentially different. 

In this paper, we investigate the phonon modes at fi- 
nite temperature in graphene and SWCNT, by using the 
mode projection technique combined with MD simula- 
tion. The carbon-carbon interaction is described by the 
second-generation Brenner inter-atomic potential. 16 Our 
results show that the phonon spectrum of flexural mode 
in graphene is renormalized from quadratic into linear 
at finite temperature. This renormalization effect should 
have significant influence in many physical phenomena, 
since the flexural mode is the first mode to be thermally 
excited due to its parabolic dispersion. The renormal- 
ization does not occur in SWCNT. This difference be- 
tween graphene and SWCNT is attributed to the two- 
fold degeneracy of the flexural mode in SWCNT result- 
ing from its quasi-one-dimensional structure. This two- 
fold degeneracy is robust against the PPS at finite tem- 
perature and it can help to protect the quadratic prop- 
erty of the phonon spectrum for flexural mode. We also 
study optical phonon modes in graphene and SWCNT. 
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The frequency of the out-of-plane optical phonon mode in 
graphene decreases linearly with the increase of temper- 
ature with a slope of 7.5 x 10 -3 cm -1 /K. The life time of 
this optical mode is shorter at higher temperature, since 
the PPS is stronger at higher temperature. For the two 
degenerate in-plane optical phonon modes in graphene, 
we find that the frequency shows increasing behavior in 
low temperature region, and will decrease with further 
increase of temperature. The frequency reaches a maxi- 
mum value around 600 K. In wider SWCNT (n, n) with 
n > 6, the axial optical mode behaviors similarly as the 
in-plane optical mode in graphene. However, in narrow 
tubes, the frequency of the axial optical mode decreases 
monotonically with the increase of temperature. The life 
time of the axial optical phonon mode is not very sen- 
sitive to the diameter, and they can be well described 
by a simple formula for SWCNT with all diameters and 
graphene as: t(T) = 560/T, which is quite comparable 
with what the experimentalists have observed. 

The rest of the paper is organized as follows. In Sec. II, 
we present the detailed calculation procedures. Sec. Ill 
is devoted to the main results and relevant discussions. 
The conclusion is in Sec. IV. 



II. CALCULATION PROCEDURE 
A. graphene and SWCNT structures 

The coordinate system in graphene is given in 
Fig. Q] (a). The two-component basis vectors of the hon- 
eycomb lattice structure are^ a\ = (V3/2,-l/2)a and 
a-2 = (V3/2,l/2)a, with lattice constant a — 2.46 A. 
Correspondingly, the two basis vectors in the reciprocal 
space shown in Fig. □ (b) are b x = (1/2,-^3/2)6 and 
b-2 = (1/2, v3/2)&, where the reciprocal lattice constant 
b = 47r/v3a- In lattice structures, there are some ar- 
bitrariness in the choice of basis vectors. The only re- 
quirement is that the area (two-dimensional lattice) or 
the volume (three-dimensional lattice) of the unit cell 
defined by the basis vectors does not changed At the 
end of this part, we will note that a pair of unsymmet- 
rical basis vectors are more convenient for us, instead of 
a\ and a,2 which have reflection symmetry with respect 
to the x axis. As shown in Fig. Q] (c), we choose ba- 
sis vectors as A% = a\ + a,2 and A2 = 0,2- The A\ and 
A2 can satisfy the requirement for basis vectors, because 
A\ x A2 — a\ x 3,2- Corresponding to A\ and A2, the two 
reciprocal basis vectors are B\ = b\ and B2 = — &i + &2 
as shown in Fig. [T] (d) . 

Fig. Q] displays that both lattice and reciprocal lattice 
of graphene sheet have the same honeycomb structure 
with D§h group symmetry^ In the reciprocal lattice 
space, high symmetry points T, M and K are impor- 
tant in many physical phenomena. Usually, the phonon 
dispersions along high symmetry lines between these 
points are shown. The Y point is in the center of the 




FIG. 1: (Color online) (a) and (c) are two sets of basis vectors 
in graphene. (b) and (d) are the corresponding reciprocal ba- 
sis vectors. A simulation sample of graphene with (n, m)=(18, 
5) is shown in (e). 

Brillouin zone with zero wave vector. The M point 
is at YM = 6i/2 = Bi/2. The K point locates at 
YK = (—61 + 62)/3 = B2/3. Another important wave 
vector is YKM = (-61 + & 2 )/2 = B 2 /2. 

A piece of graphene can be denoted by a set of two in- 
tegers (n, to). They define two lattice vectors R\ = nA\ 
and i?2 = mA^, which form two neighboring edges of the 
graphene. In this graphene sheet, there are n unit cells 
repeated in A\ direction and to unit cells in A2 direction. 
A total number of 6nm phonon modes can be found in 
this system. By applying periodic boundary conditions, 
the Bloch theorem guarantees that these phonon modes 
can be labeled bynxm wave vectors: k = k\B\ + ^2^2, 
where k\ — ji/n, j\ £ (— n/2,n/2] and k2 = fa/m, 
]2 € (—to/2, to/2]. The Bloch theorem connects the wave 
vector in the reciprocal lattice space to a system struc- 
ture in the real lattice space. This connection leads to a 
one to one mapping from the phonon modes we are going 
to investigate to the shape (also size) of the graphene on 
which we are going to perform MD simulation^ For ex- 
ample, to study a number of phonon modes with wave 
vectors equally distributed between Y and M, we should 
investigate a graphene sheet with n = (Nk — 1) * 2 and 
to = 1. In another situation, suppose that we are going to 
study a number of phonon modes with wave vectors 
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FIG. 2: (Color online) The vibrational displacement for the 
axial optical phonon mode in SWCNT (3, 3) and (10, 10). 



B. MD simulation details 

In our simulation, the second-generation Brenner 
inter-atomic potential 16 is implemented to describe the 
carbon-carbon interaction. Periodic boundary conditions 
are applied in all growth directions. The Newton equa- 
tions of motion are integrated within the fourth order 
Runge-Kutta algorithm, in which a time step of 0.2 fs 
is used. The micro-ensemble is applied. Before the MD 
simulation, the structure is optimized to the equilibrium 
position at zero Kelvin. Then the initial position and 
velocity are randomly assigned according to Gauss dis- 
tribution. 5 x 10 4 simulation steps are used for the system 
to reach thermal equilibrium state. The trajectory data 
are recorded for 5 x 10 5 simulation steps. The variation 
in the temperature is within 2%. We have also run MD 
simulation for 1 x 10 6 simulation steps, and the differ- 
ence introduced is ignorable. The largest system in our 
simulation is SWCNT (20, 20) with N = 178, which has 
14240 carbon atoms. 



C. mode projection method 



equally distributed along TKM , we have to do simulation 
on graphene with n — 1 and to = (iV/. — 1) * 2. It is appar- 
ent that larger system is in need if we want to investigate 
more phonon modes. Fig. [T] (e) shows the structure of a 
graphene sheet with (n, m)=(18, 5). This sample can be 
used to study those ten phonon modes with wave vectors 
equally distributed along TM. We have technically set 
to to be 5 instead of 1, because it should be better to 
ensure the sizes of the graphene in all directions larger 
than the interaction range of the Brenner potential. 

It is important to note that we can investigate phonon 
modes along all high symmetry lines (both TM and 
TKM) with varying integers (n, to) by using the two 
unsymmetrical basis vectors A\ and A-i . If the two sym- 
metric basis vectors a\ and 3,2 are used, the integers n and 
to are equivalent to each other; thus only those phonon 
modes distributed along the high symmetry line TM can 
be studied. 

The SWCNT (n, n) is a quasi-one-dimensional 
system.— The coordinate system in SWCNT is shown in 
Fig. [21 We have set the x direction to be the tube axis. 
SWCNT is periodic in the x direction with basis vector 
A = ae Xl and there are An carbon atoms in the unit cell. 
A general lattice vector is denoted by only one integer I 
as R{1) = IA. Correspondingly, the Brillouin zone is a 
segment B = 2tt /ae x . An arbitrary wave vector can be 
described by k = kB, with k = j/N, j € {-N/2, N/2], 
where N is the number of unit cells. Similar as the sit- 
uation in graphene, the Bloch theorem requires that if 
we want to study a number of Nk phonon modes with 
wave vectors k e [0,1/2], we need to carry out MD 
simulation for a SWCNT with the number of unit cell 
N = (N k - 1) * 2. 



For systems with pure harmonic interaction, phonons 
are free and have no interaction with each other. As a 
result, their frequency is temperature independent and 
the life time is infinite. This is an ideal situation at 
K, however the PPS is inevitable at finite temperature in 
real materials like graphene or SWCNT. The PPS leads 
to the renormalization of the frequency and gives a finite 
life time for the phonon. Since the PPS is sensitive to 
the temperature, the phonon frequency and life time will 
also have temperature dependence. This nonlinear effect 
on the phonon modes can be investigated by the mode 
projection method. 13,19 Following are the three steps to 
implement this method, and these are also the three steps 
in our numerical procedure. We demonstrate this method 
in graphene. For SWCNT, all processes are the same 
except some notation substitutions. 

(i) . We run MD simulation and record the trajectory 
profile for each atom (Ji,l2,s), say f[h,l2,s). The two 
integers l\ and I2 are lattice indexes for each unit cell in 
graphene with lattice vector as Rili^h) = I1A1 + I2A2. 
The index s = A, B refers to the two nonequivalent car- 
bon atoms in each unit cell. The relative position for s 
in the unit cell is r s . The equilibrium position of atom 
(h,h, s) at K is r (li, l 2 , s) — R(li,l2) + r s . As a result, 
the vibrational displacement can be obtained through 
u(li,l2, s) — r{l\, hi s) — ro(h, h, s). With the evolution 
of MD simulation, time dependent vibrational displace- 
ments are achieved u(l±, I2, s, t). 

(ii) . In periodic lattice structures, the Bloch theorem 
can be applied^, 

u(hi 2S ) = ^^EE4 <r)eiS ** I,a ^ ,r) (*i 00a )' (1) 



4 



where each phonon mode is denoted by two indexes (fc, 
a), k is the wave vector and a is the branch index for 
phonons at same wave vector k. N± and N2 are the num- 
ber of repeated unit cells in two basis vectors' directions. 

is the normal mode operator, and ^ a \k\00s) is the 
polarization vector for this phonon mode. Considering 
the orthogonality and completeness of the polarization 
vectors, the above expression can be reconstructed into 
another equivalent form: 

Qf = e-*^^ T) (fc|00 S )t • u(hl 2 s), (2) 

hhs 

where the vibrational displacement utlil^s) has been ob- 
tained in step (i) by MD simulation. In this step (ii), we 
get the normal mode operator for each phonon (fc, a). 

(hi). To extract the properties of each phonon mode 
from its normal mode operator, the normalized auto- 
correlation function g(t) is calculated: 

Q (t) = fc fc (3) 

<Q<>> (0)QW(0)*}' 

For free phonons, this equation simply gives g(t) — 
coswi, which oscillates without decay indicating an in- 
finite life time. The PPS will lead to finite phonon 
life time. Under the single mode relaxation time 
approximation^ - — we have g(t) = coscjie - */ 7 ", where 
lj is the renormalized frequency and r is the life time. 

Briefly speaking, the MD simulation at finite tempera- 
ture T is carried out to record the information of normal 
mode operator for each phonon (k, a). The PPS has 
been considered in MD simulation through the nonlinear 
Brenner potential. Then the correlation function g(t) is 
calculated and fitted to g(t) = coswte~*/ r . In this way, 
we obtain the phonon frequency and life time at temper- 
ature T. 



III. CALCULATION RESULTS AND 
DISCUSSION 

A. phonon modes in graphene 

Fig. [3] illustrates the process of extracting the phonon 
frequency and life time from the correlation function g{t). 
Panel (a) shows correlation functions at 300 K (blue on- 
line) and 1000 K (green online). g(t) decays faster at 
1000 K than that at 300 K, indicating shorter life time of 
this phonon mode at higher temperature. Panel (b) is the 
fitting of the correlation to a function g(t) = cos wte~*/ T . 
The fitting procedure is of high accuracy since the fitted 
curve and the MD simulation results are almost indistin- 
guishable. This means that the single mode relaxation 
time approximation is quite reasonable. In our numeri- 
cal fitting procedure, the relative uncertainty for the fre- 
quency uj (life time r) is about 10 -6 (10 -3 ). The first 




FIG. 3: (Color online) (a) is the normalized correlation func- 
tion g(t) at 300 K (blue online) and 1000 K (green online). 
The g{t) at 1000 K is fitted in (b) and (c). 



quarter of panel (b) is enlarged in (c) to highlight the 
details of the fitting. 

In graphene, the acoustic phonon spectrum in the z di- 
rection shows a quadratic dependence on wave vector k: 
uj(k) oc k 2 , which is the so-called flexural mode. Fig. 0] 
compares the phonon spectrum of flexural mode at 0, 
300 and 1000 K. The phonon spectrum at K is calcu- 
lated directly by solving the eigen value problem of the 
dynamical matrix. The phonon spectrum is renormal- 
ized obviously, and the curve shows an up-shift at room 
temperature compared with that at K. Because of this 
up-shift in the phonon spectrum, it is possible that the 
dispersion will be renormalized from quadratic into lin- 
ear. Indeed, such behavior is observed as shown in the 
inset of Fig. [4j which shows the phonon spectrum in the 
limit of long wave length. It shows clearly that the room 
temperature phonon spectrum can be fitted to a linear 
function u(k) — ck with slope c = 79, which corresponds 
to an acoustic velocity of 0.4 kms -1 after converting into 
standard units. 

This renormalization of flexural mode at finite temper- 
ature can take effect in almost all thermal phenomenon, 
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FIG. 4: (Color online) Renormalization of the flexural mode 
in graphene. Inset is the linear fitting of the phonon spectrum 
in long wave length limit. 



as it has the lowest frequency and will be firstly excited 
with increasing temperature. In some physical processes, 
the quadratic property of the phonon dispersion is im- 
portant. It is necessary to consider this renormalization 
effect at finite temperature in those physical processes. 
For example, the flexural mode has important contribu- 
tion to the resistance in low temperature region, and the 
electron scattering rate due to the flexural mode 6 will not 
diverge with temperature and need to be further renor- 
malized. In the phonon thermal conductance, after the 
renormalization of the phonon spectrum, the contribu- 
tion of flexural mode to the thermal conductance will 
depend on temperature as k oc T 2 instead of the T . 

Now we focus on the discussion of the three optical 
phonon modes at T point in graphene. Fig. [5] shows the 
frequency and life time of the optical mode in z direction 
at different temperatures. From panel (a), the frequency 
decreases monotonically with the increase of temperature 
and the slope of the decreasing curve is about 7.5 x 10~ 3 
cm _1 /K. In Panel (b), an overall decrease in the life time 
is revealed, however the variation is a little bit for this 
phonon. Fig. [6] shows the temperature dependence of 
the frequency and lifetime for the two in-plane optical 
phonon modes, which are degenerate due to the D§h 
point group in graphene. The degeneracy of these two 
optical phonon modes are well preserved, although the 
value of frequency and life time vary at different tem- 
perature. This preservation of degeneracy implies that 
the Dsh symmetry in graphene is protected even at finite 
temperature. An interesting phenomenon is that the fre- 
quency increases with increasing temperature in low tem- 
perature region, and it shows decreasing behavior after 
reaching a maximum value around 600 K. The life time 
has smaller value at higher temperature, since the PPS 
is stronger at higher temperature. Our theoretical value 
(about 1.5 ps at 300 K) is in good agreement with the ex- 
periment result, where the life time of this mode is about 
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FIG. 5: (Color online) Frequency and lifetime vs. tempera- 
ture for the out-of-plane optical phonon mode in graphene. 
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FIG. 6: (Color online) Frequency and lifetime vs. tempera- 
ture for the two degenerate in-plane optical phonon modes in 
graphene. 
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FIG. 7: (Color online) Phonon dispersion curves for the four 
acoustic phonon modes in SWCNT (3, 3). The four acoustic 
phonon modes are: doubly degenerate flexure mode, twisting 
mode, and longitudinal acoustic (LA) mode. 



1.2 ps at room temperature and is inverse proportional 
to the temperature^ The agreement between our theo- 
retical results and the experimental value discloses two 
facts: (i) the PPS dominates the phonon life time under 
the experimental conditions which is also confirmed by 
the inverse proportional temperature dependence of life 
time, and (ii) the mechanical exfoliated graphene sam- 
ples are of high quality thus other scattering mechanisms 
are ignorable. 



B. phonon modes in SWCNT 

In this section, we investigate the phonon modes in 
SWCNT at finite temperature. The flexural modes in 
SWCNT are the two transverse acoustic phonon modes, 
which shows quadratic spectrum at K. They are de- 
generate because of the quasi-one-dimensional structure 
where the two vertical directions (y and z) are equiva- 
lent to each other. Fig. [7] illustrates the four acoustic 
phonon modes in SWCNT (3, 3): two-fold degenerate 
flexure mode, twisting mode, and the longitudinal acous- 
tic mode. We focus on the doubly degenerate flexure 
mode. Fig. [3] shows the flexural modes in long wave 
length limit of different SWCNT (n, n) at room tem- 
perature. The phonon dispersion has larger frequency in 
wider tubes. It is quite interesting that the quadratic 
property of the phonon spectrum for all tubes studied 
is not broken by PPS at finite temperature, which con- 
firms previous results^ The inset of the figure compares 
the phonon spectrum for the flexural mode in SWCNT 
(3, 3) at 0, 300 and 1000 K in a larger range of wave 
vector. It shows that the phonon spectrum has a down- 
shift at room temperature compared with K. This fact 
guarantees that the phonon spectrum will not be renor- 
malized into linear dispersion. This result is opposite 



FIG. 8: (Color online) The flexural mode in SWCNT (n, n) 
at room temperature. From bottom to top, n=3, 6, 10, 12, 
14, 16, 18 and 20. Inset shows the overall figure of the flexural 
mode in tube (3, 3) at 0, 300, and 1000 K. 



to what we have observed above in graphene, where the 
flexural mode is renormalized into linear dispersion at 
finite temperature. One possible reason for this differ- 
ence between graphene and SWCNT is that the flexu- 
ral mode in SWCNT is two-fold degenerate while it is 
non-degenerate in graphene. The degeneracy in phonon 
spectrum is robust against PPS at finite temperature, 
which is also the case for the two degenerate in-plane 
optical phonon modes in graphene as we have discussed 
in previous section. This degeneracy may have impor- 
tant effect in the protection of the quadratic property 
of the phonon spectrum of flexural mode. The degen- 
eracy originates in the quasi-one-dimensional structure 
of the SWCNT, where a dynamical band gap is opened 
for the flexural mode due to the four-phonon process.— 
This dynamical gap may facilitate the protection of the 
quadratic phonon spectrum of flexural mode in SWCNT. 

Fig. [5] exhibits the vibrational morphology of the opti- 
cal phonon modes in the tube axial direction in SWCNT 
(3, 3) and (10, 10). In large diameter limit, this mode 
turns to be the in-plane optical phonon mode in the 
graphene. The frequency of this mode is shown in 
Fig. 12] (a), where the in-plane optical mode in graphene 
is also shown for comparison. The frequency is smaller in 
wider SWCNT and approaches to the value of graphene 
in the limitation of large diameter. In wider SWCNT 
(n, n) with n > 6, the curve behaviors similarly as 
that shown above for graphene. The frequency increases 
with the increase of temperature at low temperatures, 
and then decreases with further increasing temperature 
above 600 K. The maximum value of frequency is reached 
around 600 K in all wide tubes. Situation changes in 
narrow SWCNT where the frequency decreases mono- 
tonically with the increase of temperature. Fig. [9] (b) 
shows that the life time of this mode is not very sensitive 
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FIG. 9: (Color online) Frequency and life time v.s tempera- 
ture for the axial optical phonon mode in SWCNT. 



to the diameter, and it is roughly inversely proportional 
to the temperature. All data of life time can be fitted 
to a function r(T) = 560/T. The experimental value of 
the life time for this optical phonon at room temperature 
is around 1.2 ps, varying in samples prepared by differ- 
ent methods i 22 1 23 The measured life time is insensitive to 
the diameter or chiral property of the SWCNT, and it 
is inverse proportional to the temperature in the exper- 
imental range [300, 450]K. The experimental value (1.2 
ps) for the life time is smaller than our theoretical results 
(1.8 ps). This difference is because of other inevitable 



scattering mechanisms, such as defects or impurities, in 
the experimental SWCNT samples. Actually, these scat- 
tering mechanisms have exposed themselves through the 
variation in the value of life time for samples prepared 
by different methods. 



IV. CONCLUSION 

To conclude, we study the frequency and life time of 
different phonon modes at finite temperature in graphcne 
and SWCNT, by using the mode projection technique 
combined with MD simulation. We pay special attention 
to the renormalization of the quadratic phonon spectrum 
of the fiexural mode. In graphene, the phonon spec- 
trum of fiexural mode is renormalizcd into linear at finite 
temperature. The importance of the renormalization is 
illustrated by some examples. However, the renormal- 
ization does not happen in the SWCNT, because of the 
two-fold degeneracy of the fiexural mode in the quasi- 
one-dimensional SWCNT. For the three optical phonon 
modes at T point in graphene, we study the tempera- 
ture dependence for their frequency and life time. The 
frequency of the out-of-plane optical phonon mode in 
graphene decreases with the increase of temperature with 
a decreasing slope of 7.5 x 10~ 3 cm _1 /K. The two in- 
plane optical phonon modes in graphene remains degen- 
erate at finite temperature. With the increase of temper- 
ature, the frequency will increase to a maximum value 
around 600 K and then decrease. A similar phenomenon 
is found for the axial optical mode in wider SWCNT (n, 
n) with n > 6. In narrow tubes, the frequency of the 
axial optical mode decreases monotonically with the in- 
crease of temperature. The life time of the axial optical 
phonon mode in all SWCNT studied and graphene can 
be fitted to a function: t(T) — 560/T. This result is in 
good agreement with the experiments. 
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